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Abstract 

An accurate numerical consideration of ID spinless fermion model with 
next-nearest neighbour (NNN) interactions is carried out for the electron 
concentrations 4/7. It is shown that depending on the parameters of 
the model it can be either Luttinger liquid or bipolaron liquid. In the 
former case competing interactions can result in a smooth behavior of 
one-electron distribution function at the Fermi surface with a divergence 
in the second derivative with respect to the quasimomentum (and not 
in the first one, as usual). In this connection, uncommon photoemission 
spectra for some ID conductors can be explained qualitatively. 

PACS 71.10.Pm, 71.10.Fd 

Frustrations and competing interactions are very important for strongly cor- 
related electron systems. As it was shown in a seminal work|l| they can change 
drastically a ground state of quantum antiferromagnets leading to the forma- 
tion of "resonating valence bond" (RVB), or quantum spin liquid, state (see 
also Q and references therein). This state has been found later in a series of 
compounds with both magnetic (spin liquid) and charge (pseudospin liquid) 
degrees of freedom, see, e.g., the discussion in Refs.|^, Q].The effects of com- 
peting interactions were investigated recently for one-dimensional (ID) spinless 
fermion model p, 0, u\. Usually this model is used for the description of Verwey 
metal-insulator transitions and charge ordering phenomena in Fe304, Ti407 
and other d- metal compounds M M 0|. It is expected that in the half- filled case 
the growth of the Coulomb interaction leads to the transition from the metal- 
lic state to the insulating charge-ordered one. It appeared, however, that the 
competing interactions can result in the stabilization of metallic phase without 
charge ordering for arbitrarily strong Coulomb inter action fi 0], this phase being 
not usual (for ID systems) Luttinger liquidH. Here we investigate this model 
beyond half-filling. 

We proceed with the Hamiltonian 

L L L 

H ^ -t^^ \c\ci+i -I- 4+iCiJ + V^ ^ n^n^+i + V'^^ nin^+2 (1) 

Z— 1 2—1 Z— 1 
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where cl,Ci are Fermi creation and annihilation operators on site i, Ui — c\ci. 
Further we use periodic boundary conditions (L is the length of the ring) and 
put i = 1 at the presentation of the computational results. 

Earlier we investigated the cases p — 1/2 and p = 2/3 where p = N/L is 
the electron concentration, TV is the number of electrons^, M. In that cases the 
ground state is charge-ordered and insulating for small enough t, except special 
parameter ratios {V = 2V' for p — 1/2 and y = or y = for p — 2/3) 
when the ground state has macroscopically large degeneracy g = exp(LS'), S is 
the entropy per cite in the limit t — 0. Consider now a case 1/2 < p < 2/3. 
It is easy to demonstrate that S is finite at i = and arbitrary positive V 
and v. Therefore one can expect that the ground state is gapless (conducting) 
at any small but non-zero t. However, the character of the electron motion 
turns out to be different for V > 2V' and V < 2V' . In the former case the 
processes of the motion of a single electron without the change of the energy 
of Coulomb repulsion are possible but in the latter one the electrons can move 
only as a coupled pairs on neighbouring sites (bipolarons). For example, the 
configurations [...011011001...] and [...011001101...] which are distinguished 
by the motion of one bipolaron in right direction have the same interaction 
energy (here 1 and are labels for occupied and empty sites). At the same 
time, the configurations [...010111001...] or [...011010101...] have the additional 
energy equal to V and 2V' — V, correspondingly. The latter is the value of 
bipolaron binding energy at t = 0,V' < V < 2V' . Note that this bipolaron 
state is similar to Shubin-Vonsovsky "maximum polarity" statc[0]. Thus, the 
cases V > 2V' and V < 2V' are physically diff'erent and has to be considered 
separately. We investigate both of them numerically for the case p = 4/7. 
Parameters of the ground state, correlation functions, etc. are calculated by 
the Lanczos method with the extrapolation to the infinite length of the chain 
L ^ oo] the computational methods are described in Refs.|^, ^. 

Consider first the region V > 2V' . One can assume that the system is 
Luttinger liquid (LL)||l^ in this case. There are three branches of low-energy 
excitations in LL: density fluctuation, current, and charge excitations with the 
velocities vs,vj and vn, correspondingly. The first one is connected with the 
variation of the total energy E of the system under the variation of the total 
momentum P, vs = SE/SP] the second one (wj), with the variation of the 
energy under a shift of all the particles in momenta space, and the third one, 
with the variation of the chemical potential p under the change of the total 
number of particles N, vn = (L/tt) 6p/6N. The velocity vj is connected with 
the conducting properties of the system[|3| Q: vj — 2TrD where D is the Drude 
weight (a spectral weight of the zero-frequency peak in the conductivity) . Thus 
vj 7^ is one of the criterion of the metallic state [O. As it was shown by 
Haldane|n3], the identity 

{vsf = VjVn (2) 

takes place. Therefore it is possible to introduce the parameter 

e"^*^ = vn/vs = vs/vj (3) 



which determines the time and space asymptotics of the correlation functions 
of the system. One-electron distribution function n{k) as a function of the 
quasimomentum k has a power-low singularity near the Fermi surface, k = kp 



n{k) « n(kp) — Csgn(k — kp)\k — kpl"^ 



(4) 



with a ~ cosh2</? — 1. One-particle local spectral density p (E) has a singularity 
near the Fermi energy /i with the same exponent, p (E) oc |-B — /i| which can 
be measured therefore by photocmission spectroscopy. There are various experi- 
mental data for inorganic materials with charge-density- waves [|5[ |lj, |l^ as well 
as for the organic superconductor (TMTSF)2PF6|lq] which confirms power-law 
singularity in p (E) predicted by LL theory. It is interesting that in the latter 
case the value a ~ 1.25 has been found. At the same time, in all the known 
models with short-range interactions we always have a < 1 (in particular, in 
the ID Hubbard model one can prove that a < l/8|l9|). Our numerical calcu- 
lations show that a < 1 is also in the model under consideration for one-half - 
or two-third- fillings W. Here we demonstrate that a can be larger that 1 for 
1/2 < p < 2/3 which give a qualitative explanation of the experimental data|jlq| . 
The computational results are presented in the Table. 

Table 1: Velocities vj,vn, and vs extrapolated to L — > cxd, the ratio x = 
vjVn/ (vs) , and the parameter a from Eq.(E|) for different V,V' [t ~ 1). For 
the lower part of the Table there is no Luttingcr liquid state (%, wjv — > oo at 
L ^ oo and a is not defined). 
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The finiteness of vj and, consequently, of the Drude weight proves the metal- 
lic nature of the ground state. One can see that for V > 2V' Eq.(0) is valid (but 



not for V < 2V' where the system is not LL but bipolaron liquid, see below). 
The exponent a appears to be maximum along the line V = (11/40)1^. In 
particular, for V — 120, V = 33 one has a ~ 1.68 which is even larger than the 
experimental value for (TMTSF)2PF6@. 

Consider now the case V < 2V' . It is shown from the Table that Eq.(0) is 
violated in this case, at least, for large enough V and V . As it was discussed 
above the reason is the formation of bipolarons. There are two evidences of the 
coupling of electrons in the pairs for large enough 2V' — V. First of all, the 
energy gap defined ms A = E {N + 1) + E {N - 1) - 2E (N) is finite (Fig. 1); 



Figure 1: Energy gap as a function of V' at V ^ 20 for L =14(o), 21(d), 28(a), 
and the results of the extrapolation to L ^ oo(»). 

at the same time, the Drude weight is also non-zero which means that A is the 
superconducting, or bipolaron, gap but not insulating one. Then, the ground 
state energy appears to be oscillate as a function of magnetic flux in the ring 



$ with the period corresponding to the charge 2e of current carriers. The 
corresponding results are presented in Fig. 2, 



Figure 2: Flux dependence of the total energy of the ring (L = 21) for T^ = 20 
and different V: V = 12 (diminished by a factor 2)(o); V'=U{a); y'=15(A); 
r=20(v). 

the flux being introduced by the replacement t — *■ iexp(27ri$) in Eq.(|l|). 
According to these criteria, bipolarons do not form at small enough 2V' — V^ 
t. To understand better this transition region we fit the total energy by the 
expression 

E{^) -E{0) =/3(l-cos27r$)+7(l-cos47r$) (5) 

where the first term corresponds to the "normal" current carriers and the second 
one to the bipolarons. The fitted parameters f3 and 7 are shown in Fig. 3. 

One can see that there is a coexistence of bipolaron conductivity and the 
usual one in the transition region V ~ 2V' 3> t. It would be very interesting to 



Figure 3: Fitted parameters 13/ ^J (3"^ + 7^ (•) and 7/\//3^ + 7^ (o) from Eq.(||) 
as a function of V' for V = 2Q and L = 21. 



investigate an explicit nature of this coexistence (phase separation or a homo- 
geneous state) but for this aim computations for larger clusters are necessary. 

In the limit of small t one can simplify the problem using the operator 
perturbation theory developed first by BogoliubovpO]. In the lowest order in 
t/ {2V' — V) one can obtain 

P ^ / \ 

^^// " ~ 2V'-V '^ ^ (4c»+lcLlCi + 4c^^l4+lC^) 'P ■ (6) 

i—l 

where V is the projection operator on the state with the lowest Coulomb energy. 
This Hamiltonian describes the bipolaron motion with the effective transfer 
integral t"^/ {2V' — V). Our numerical results give for the Drude weight 

vj - 27rD ^ '^ (7) 

2V' -V ^ ^ 

For small enough t the latter can be arbitrarily small. This differ drastically 
the case under consideration from the usual Verwey transition. One can show 
for p = 1/2, V' ~ Q that near the transition point V — 2t the Drude weight has 
a finite jumpQ 

vj{V ^2t-0)^7r/2; vjiV ^2t + 0) = 0. (8) 

To conclude, we investigate numerically the ground state of ID spinless 
fermion model with competing interactions beyond half-filling. Depending on 
the ratio of NN and NNN interactions it is either Luttinger liquid (with un- 
usually large exponent a which can be even larger than 1) or bipolaron liquid. 
These results may be important to understand properties of low-dimensional 
conductors; in particular, the value a > 1 has been observed in photoemission 
experiments for (TMTSF)2PF6||l|l. 
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